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Abstract 

We calculate threshold soft-gluon corrections to total cross sections and transverse mo- 
mentum distributions for bottom and charm quark production in fixed-target experiments, 
considering both pp and 7r~p interactions. We investigate the quality of the near thresh- 
old soft-gluon approximation at next-to-leading order (NLO) and calculate next-to-next- 
to-leading order (NNLO) corrections through next-to-next-to-next-to-leading-logarithmic 
(NNNLL) accuracy, including some virtual terms. We find that the NNLO threshold 
corrections reduce the factorization and renormalization scale dependence of the cross 
sections. 



1 Introduction 



The latest calculations for heavy quark hadroproduction have included next-to-next-to-leading- 
order (NNLO) soft-gluon corrections to the double differential cross section PQ |2l El IH E] from 
threshold resummation techniques [HIIZIIH]- These resummations are a consequence of the factor- 
ization properties of QCD which separate cross sections into universal nonperturbative parton 
densities and a perturbatively-calculable partonic cross section. Near threshold there is limited 
phase space for the emission of real gluons so that soft-gluon corrections dominate the cross 
section. These Sudakov corrections have the form of logarithmic plus distributions, singular at 
partonic threshold. Threshold resummation techniques organize these singular distributions to 
all orders, thus extending the reach of QCD into the near-threshold region. 

Calculations of bottom and especially charm production are still not under solid theoretical 
control. A good understanding of the bottom cross section is important for HERA-B, where b- 
quark production is near threshold. The charm cross section is of particular interest for heavy 
ion physics. Although many future heavy ion experiments will be at high collider energies, 
y/S > 130 GeV, some of the current and future experiments, like those at the SPS, are in the 
near-threshold region. The NA60 experiment will take heavy ion data at \^ =17.3 GeV and 
pA data at = 29.1 GeV. A new facility is being built at the GSI ^ that will measure charm 
near threshold with = 6.98 GeV. 

Because the charm quark mass is a few times Aqcd, it is generally treated as a heavy quark 
in perturbative QCD calculations. However, its relative lightness results in a rather strong 
mass and scale dependence for the total cross section. There is also a rather broad spread in 
the charm production cross section data at fixed target energies. 

In this paper, we increase the accuracy of previous soft-gluon calculations for bottom 
and charm production. The soft corrections that we calculate take the form of logarithms, 
[ln^(xth)/a;th] + 7 with I < 2n — 1 for the order a" corrections, where Xth is a kinematical vari- 
able that measures distance from threshold and goes to zero at threshold. NNLO calculations 
[n = 2) for bottom and charm quark production have so far been done through next-to-next- 
to-leading-logarithmic (NNLL) accuracy, i.e. for the scale-independent terms, including leading 
logarithms (LL) with / = 3, next-to-leading logarithms (NLL) with I = 2, and NNLL with / = 1 
PQ 121 13 HI- In Refs. 12 El El 5 heavy quark cross sections were studied in both single-particle- 
inclusive (IPI) and pair-invariant-mass (PIM) kinematics. Important differences between the 
two kinematics choices were found in both the parton-level and hadron-level cross sections, even 
near threshold. Thus subleading, i.e. beyond NNLL, contributions can still have an impact 
on the cross section. Their inclusion is clearly needed to bring the calculation under further 
theoretical control. The subleading terms indeed minimize the kinematics dependence of the 
top quark production cross section j^l. 

In Refs. [H 12] we studied top and bottom quark production at NNLO-NNLL. In Ref. jll 
we studied charm quark production at NNLO-NNLL. More recently we extended our calcu- 
lations for top quark production [S], using the methods and results of Ref. [TU], to include 
additional subleading NNLO soft corrections, including next-to-next-to-next-to-leading loga- 
rithms (NNNLL), with n = 2, / = 0, as well as some virtual 6{xth) corrections. We showed in 
Ref. [Sj that the subleading corrections do indeed bring the IPI and PIM results into agreement 
near threshold for both the qq — > QQ and the gg QQ channels while the discrepancies away 
from threshold are also diminished, especially in the gg channel. Thus the threshold region is 
brought under better theoretical control. In this paper, we apply these new terms to bottom 
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and charm quark production. 

In the following section, we briefly discuss the differences between the IPI and PIM kine- 
matics choices and then investigate which choice is a better approximation to the full NLO 
result. We flnd that IPI kinematics is a far better choice. Thus we only calculate the cross 
section in this kinematics in the remainder of the paper. 

Section 3 discusses bottom quark production in 7r~p and pp collisions, particularly at HERA- 
B. We present the hadronic total cross sections for a range of energies near threshold. We also 
calculate the transverse momentum distributions for bottom production at HERA-B. In Section 
4 we study charm production in 7r~p and pp collisions. Our studies focus on the kinematics 
of the proposed GSI facility and the CERN SPS proton and ion fixed target programs. We 
conclude with a summary in Section 5. 

2 Singular distributions and kinematics dependence 

As we discussed in the introduction, the soft-gluon corrections that we calculate take the form 
of logarithms, [ln^(xth)/a;th] + , where kinematical variable that measures the distance 

from partonic threshold. The exact definition of x^^ depends on the kinematics we use for the 
calculation of the cross section. 

We study the partonic process ij — *• QQ where Q is the produced heavy quark and ij can 
be either qq or gg. A more detailed discussion of the kinematics can be found in Ref. |2J. 

In IPI kinematics, a single heavy quark is identified so that 

^{Pa)+3{p,) Q{pi)+X\Q]{P2) (2.1) 

where Q is the identified bottom or charm quark of mass m and X[Q] is the remaining final 
state that contains the heavy antiquark Q. We define the kinematical invariants s = {pa+PbY, 
ti = ipb — PiY "~ ""^^5 ""i = iPa — PiY ~ ""^^ S4 = s + ti + Ml. At threshold, S4 — and 
the soft corrections appear as [In'(s4/m^)/s4] + . We note that the virtual corrections multiply 

5(S4). 

In PIM kinematics, we have instead 

^{Pa)+3{Pb) QQip)+X{k) . (2.2) 

At partonic threshold, s = M^, is the pair mass squared, ti = — (M^/2)(l — PmcosO) and 

Ui = — (M^/2)(l + /5m cos 6*) where Pm = 1^1 — Am? /M"^ and 6 is the scattering angle in the 

parton-parton center-of-mass frame. The soft corrections appear as [ln'(l - z)/{l - z)]+ with 
z = M'^/s — 1 at threshold. We note that the virtual corrections multiply 6{1 — z). 

The effect of the soft-gluon corrections varies depending on the kinematics choice because of 
uncalculated subleading terms. Thus the quality of the soft-gluon approximation, the extent to 
which the soft-gluon corrections approximate the full corrections at a given order in a.^, depends 
on the kinematics choice even though the full corrections do not depend on the kinematics. 

We thus first test the soft-gluon approximation in IPI and PIM kinematics at NLO. In 
Fig. [T]we present the ratios of the exact NLO cross sections, ctnlo, to the approximate soft-plus 
virtual NLO cross section, o"^lo' (the Born result plus the (!?(tt^) soft and virtual corrections) for 
bottom quark production in pp collisions over a range o^^/S near bottom production threshold. 
In our calculations we use the MRST2002 NNLO [TT] parton densities. We show the ratios 
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in both IPI and PIM kinematics for each of the two partonic channels separately as well as 
for their sum for three scale choices: fi = m/2,m, and 2m. Note that the factorization and 
renormalization scales are set equal and are both denoted by /i. 

It is clear from Fig.^that the soft-plus- virtual result is a much better approximation to the 
exact NLO total cross section in IPI kinematics than in PIM kinematics over the entire range 
of energies and scales shown. The same conclusion is reached for charm production, shown 
in Fig. 121 for pp collisions with fi = m and 2m. Since the IPI kinematics choice is a better 
approximation to the exact NLO cross section, here we choose to work in this kinematics only. 

One-particle inclusive kinematics is a better approximation to the total bb and cc production 
cross sections because the gg channel is dominant. For fi > m, PIM kinematics gives a better 
approximation to the NLO qq cross section, not unexpected due to the simple color structure of 
this s-channel process. However, with its more complex color structure, the gg channel is more 
amenable to the IPI kinematics choice. We also found that, when the gg channel dominates 
production, the IPI scale dependence remains small while the PIM kinematics choice has a 
large scale dependence, larger than the exact NLO [SI. 
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Figure 1: The ratios of the NLO exact, (Tnlo over the NLO soft-plus-virtual, ct^lq) cross 
sections for bottom quark production with m = 4.75 GeV are shown for the gg (dashed) and 
qq (dot-dashed) channels separately, along with their sum (solid). 

Finally, we investigate the validity of the threshold approximation if only the soft terms 
are included and the virtual ones are not. Figure El compares the ratios, 0"nlo/<7nlO) with 
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Figure 2: The ratios of the NLO exact, ctnlo over the NLO soft-plus-virtual, ct^lo' cross 
sections for charm quark production with m = 1.5 GeV are shown for the gg (dashed) and qq 
(dot-dashed) channels separately, along with their sum (solid). 



the soft gluon terms alone, to ctnlo/o'nlo bottom quark production in IPI kinematics with 
fi = m = 4.75 GeV. We see that even though the ratios with the NLO soft gluon approximation 
alone are somewhat farther from one, they are still a rather good estimate of the exact NLO 
result. 

In the following sections we present cross sections for bottom and charm quark production. 
The LO and NLO results are exact. The NNLO results are approximate and include soft- 
gluon contributions only. At NNLO we give results both for NNLL accuracy, i.e. for the 
scale-independent terms, including all In'^(s4/m^)/s4, In^(s4/m^)/s4, and In(s4/m^)/s4 terms, 
and for NNNLL+(^ accuracy, which includes, in addition to the NNLL terms, all I/S4 terms 
(NNNLL) as well as some ( terms in the virtual corrections. For full details, see Ref. [5". We 
note that the NNLO-NNNLL+(^ result includes all NNLO soft and virtual terms proportional to 
the factorization and renormalization scales, while the NNLO-NNLL calculation includes in the 
virtual contribution only scale-dependent terms proportional to the squares of scale-dependent 
logarithms. 



3 Bottom quark production 

The main difficulty in comparing the calculated QQ cross sections to data lies in the fact that 
charm and bottom quarks hadronize before decaying. Top quark production is much cleaner 
in this respect since the top quark decays before it can hadronize and the measured cross 
sections can be compared directly to theory. No attempt has been made to correct B meson 
measurements to b quark cross sections here. 
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Figure 3: The ratio of ctnlo/c"nlo; "with the soft gluon terms alone, and ctnlo/cnlo; including 
the virtual terms, are shown for bottom quark production at /i = m = 4.75 GeV. We show the 
results for the gg channel (long-dashed) with (lower) and without the virtual terms (upper); 
the qq channel (dashed) with (upper) and without the virtual terms (lower); and their sum 
(solid) with (lower) and without the virtual terms (upper). 
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3.1 6-quark production in tt j9 interactions 

There is not much data on bottom quark hadroproduction at fixed-target energies. Some of 
the earhest data were from Tr~A interactions where A is a nucleus. A hnear A dependence 
was typically assumed. The data [121 [13 [HI [HI [IH] are shown in Fig. [D along with our 
calculations. The NAlO point at = 19.7 GeV is a compilation of data taken at 140, 194 
and 286 GeV on tungsten targets measuring trimuons J^. They quote a BB cross section 
essentially independent of energy. The xp distribution of the B mesons was assumed to be 
uniform around xp = with a cross section proportional to (1 — Assumptions regarding 

the shape of the distribution and the decay branching ratios almost certainly overestimate the 
cross section. The other data is analyzed using the B decay to J/ip. The WA78 data, taken 
with 320 GeV vr" beams on a uranium target, used several hypotheses regarding production 
models and B^ — B mixing angles. The data point shown assumes the largest {xp) and the 
maximal mixing. It is the smallest cross section reported in Ref. JH]- The WA92 data were 
taken at 350 GeV on a copper target Jl] . The E653 data was taken with a 600 GeV beam on 
emulsion while the E672 data at 515 GeV was taken on a beryllium target ^E]- Only this 
latter point used a NLO calculation to extrapolate from Xi? > 0, where the data was taken, 
to all Xp. All these data rely on relatively few events and some of the assumptions should be 
taken with a grain of salt. 

Figure [D presents the energy dependence of the production cross sections as functions of 
both bottom quark mass and scale. We use the GRV98 HO proton parton densities [Ej with 
the GRS pion densities [IHI. The GRS pion densities, produced in 1999, are the most recent 
set of pion densities and are compatible with the GRV98 densities. We do not use other proton 
densities because there are no recent equivalent pion sets. 

Since the bottom quark mass is relatively large, we vary the scale between m/2 and 2m 
in Fig. [U For our central value of the bottom quark mass, m = 4.75 GeV, we present the 
exact NLO cross section (solid curve), the IPI NNLO-NNLL cross section (dot-dashed) and 
the IPI NNLO-NNNLL+C cross section (dashed). We also show the IPI NNLO-NNNLL+C 
cross sections for m = 4.5 GeV (dotted) and 5 GeV (dot-dot-dot-dashed). Because we focus 
on the NNLO-NNNLL-FC results here, we do not present the NLO and NNLO-NNLL cross 
sections for the lower and upper mass values in Fig. 0] 

It is difficult to quantitatively compare the calculations to the data because the data are 
not really compatible with each other. Therefore our discussion is only on a qualitative level. 
Note that, at m = 4.75 GeV, the IPI NNLO-NNNLL-^C result is intermediate to the exact 
NLO and the IPI NNLO-NNLL cross sections. The subleading terms reduce the overall NNLO 
corrections, resulting in a value near the average of the IPI and PIM NNLO-NNLL cross 
sections. The m = 4.5 GeV IPI NNLO-NNNLL-t-^ cross section is typically equivalent to 
[fi = m/2) or larger than [fi > m) the m = 4.75 GeV IPI NNLO-NNLL cross section while 
the 5 GeV result is similar to the exact NLO cross section for m = 4.75 GeV. Changing the 
bottom quark mass between 4.5 and 5 GeV changes the IPI NNLO-NNNLL-|-(^ cross sections 
by a factor of two. 

Near threshold, 7r~p bb production is dominated by the qq channel, due to the valence 
quark- valence antiquark contribution at large momentum fractions jTHj. The qq channel gives 
the largest contribution to the bb cross section for < 40 GeV, the range shown in Fig. [H 

We now discuss the convergence properties of the subleading terms. Various definitions of 
the 'first-order K factor', K''^^ = 0"nlo/<^lO) were discussed in Ref. [20]. In Ref. [3 [1], we 
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Figure 4: The energy dependence of bb production in vr p collisions with (a) fi = m/2, (b) 
/i = m and (c) fi = 2m, calculated with the GRV98 HO proton densities and the GRS pion 
densities. We show the NLO (sohd) and IPI NNLO-NNLL (dot-dashed) resuhs at m = 4.75 
GeV. The IPI NNLO-NNNLL+C results arc shown for m = 4.75 (dashed), 4.5 (dotted) and 5 
(dot-dot-dot-dashed) GeV. 
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Figure 5: The /^-factors for bb production in TT~p collisions with m = 4.75 GeV. We present 
Kq (solid), K^'^^ (dashed) and (dot-dashed) results for (a) n = m/2, (h) jj, = m and (c) 
II = 2m. 
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compared two of these definitions, Kq , where (Tlo "was calculated with NLO parton densities 
and a two- loop evaluation of ag, and K2'\ where ctlo was calculated with LO parton densities 
and a one-loop evaluation of a^. In both cases, the Born and 0{a^) contributions to ctnlo 
were calculated with NLO parton densities and two-loop evaluations of a^. The first definition, 
Kq'^\ indicates the convergence of terms in a fixed-order calculation while the second, K2^\ 
indicates the convergence of the hadronic calculation towards a result. Since a NNLO set of 
parton densities has recently become available, in this paper, we will use both NLO and NNLO 
parton densities to calculate Kq^\ Thus when the NNLO densities are used, ctlo in -^o^'' is 
calculated with the NNLO parton densities and a three-loop evaluation of as- The definition 
of Kq^^ then indicates the convergence of terms in a fixed-order calculation for a given set of 
parton densities. We compare Kq^^ to the IPI NNLO K factors, i^*^^^ = ctnnlo-nnll/ctnlo and 
-^sub ~ o'nnlo-nnnll+c/^nlo; the latter includes the new subleading terms. In each case, the 
NNLO K factors are either calculated with NLO or NNLO parton densities and a two-loop or 
three-loop evaluation of as at each order. If K^"^^ or K^'^^ is less than Kq'\ then convergence of 
the perturbative expansion is indicated. 

The results for 7r~p —>■ bb production are shown in Fig. El as a function of energy. We present 
the K factors for m = 4.75 GeV only since the results are very similar for the other masses. 
All the K factors are calculated with the GRV98 HO proton parton densities and the GRS 
pion parton densities. Note that i^Q^^ has the strongest energy dependence with a minimum 
around = 20 GeV, higher at lower and slowly increasing with ^/S above 20 GeV. We 
find K^^^ ^ 1 for /i = m/2, increasing to ^ 2.1 at /i = 2m. The IPI NNLO-NNLL K factors, 
K^^\ the dashed curves, are flatter but K'^'^^ is larger than Kq^^ when /i = m/2. However, K^'^'' 
is nearly independent of scale and energy. Including the subleading terms increases the scale 
dependence of i^g^^ slightly but K^'^^ < K^'^^ for all energies, K^^^ ^ 1.1 for /i = m/2 and ^ 1.4 
for = 2m. Note that i^^^-' and are both nearly independent of except for < 20 
GeV. The small qq NNLO corrections reduce the increase of the K factor at low v^- 

Another measure of convergence is the scale dependence. Going to higher orders should 
make the result increasingly independent of scale. We found this to be true for tt production 
at the Tevatron 0. Here we show the scale dependence as a function of energy for the ratio 
o"(/i = m)/cr(/x = hq) where = m/2 in the upper plot and 2m in the lower. Since cr(/i = 
m) < a{fi = m/2), the upper ratio is less than unity. Thus an increasing independence of scale 
would make the ratio approach unity from below as the cross section is calculated to higher 
orders. The scale dependence is indeed decreasing for the IPI NNLO-NNNLL+(^ cross section 
over all energies. This is a significant improvement over the NNLO-NNLL result alone, further 
from unity than the NLO ratio. Improvement is also seen for /io = 2m where the ratio should 
approach unity from above. 

3.2 6-quark production in pp interactions 

We now turn to pp production of bottom quarks. So far, three experiments have reported 
the bb total cross section in proton-induced fixed-target interactions |221I2SII23) similar 
energies. The cross sections are all based on relatively small event samples. There are two 
measurements at 800 GeV, the E789 result, using a gold target j22], and the E771 result, using 
a silicon target (23 • There is a rather large difference between the two reported cross sections. 
The E771 measurement is in agreement with the HERA-B measurement at 920 GeV, taken on 
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Figure 6: The scale dependence of bb production in tt p collisions with m = 4.75 GeV. We give 
the ratios a{fx = m)/a{n = no) for the LO (dotted), NLO (sohd), IPI NNLO-NNLL (dashed) 
and IPI NNLO-NNNLL+C (dot-dashed) cross sections. The results with jiQ — m/2 are given 
in (a) while the results with //q = 2m are given in (b). 



11 



carbon and titanium targets [25 • All three experiments use the J/ip decay channel. 

The three data points are compared to our calculations with the GRV98 HO proton parton 
densities in Fig. 13 and with the MRST2002 NNLO parton densities in Fig. |S] The results are 
again shown for the exact NLO, IPI NNLO-NNLL and IPI NNLO-NNNLL+C cross sections 
at m = 4.75 GeV as well as the IPI NNLO-NNNLL+C cross sections at m = 4.5 and 5 GeV. 
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Figure 7: The energy dependence of hh production in pp collisions with (a) fi = m/2, (b) /i = m 
and (c) /i = 2m, calculated with the GRV98 HO proton densities. We show the NLO (solid) and 
IPI NNLO-NNLL (dot-dashed) results at m = 4.75 GeV. The IPI NNLO-NNNLL+C results 
are shown for m = 4.75 (dashed), 4.5 (dotted) and 5 (dot-dot-dot-dashed) GeV. 

The gg channel dominates bb production in pp collisions over the entire energy region shown. 
Note that the pp —>■ bb cross section is considerably smaller than the n~p — > bb cross section in 
the near threshold region. The n^p cross section is larger because of the high valence quark- 
valence antiquark luminosity, absent in pp collisions. At higher energies, far above threshold, 
the cross sections in both processes become more similar. Due to the difference in the pion and 
proton gluon distributions, however, they do not become equal. 

The trends in the calculated cross sections are similar to those of Fig. HI Typically, the IPI 
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Figure 8: The energy dependence of bb production in pp collisions with (a) /i — m/2, (b) 
H = m and (c) /i = 2m, calculated with the MRST2002 NNLO proton densities. We show the 
NLO (solid) and IPl NNLO-NNLL (dot-dashed) results at m = 4.75 GeV. The IPl NNLO- 
NNNLL+(' results are shown for m = 4.75 (dashed), 4.5 (dotted) and 5 (dot-dot-dot-dashed) 
GeV. 
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NNLO results with the MRST2002 NNLO parton densities are somewhat larger than those 
with the GRV98 HO parton densities although the exact NLO results are quite similar. The 
difference is due to the somewhat larger value of A4 for the MRST2002 NNLO densities. The 
HERA-B collaboration [211 reported that with m = 4.75 GeV, fi = m and the CTEQ5M 
densities [2Ij, the average of the IPI and PIM NNLO-NNLL cross sections in Ref. |2] was 
in good agreement with their data. Since A4 for the CTEQ5M densities is larger than the 
corresponding values for either the GRV98 or MRST2002 densities, the NNLO contribution is 
larger for the CTEQ5M densities. The IPI NNLO-NNNLL+C result, similar to the IPI and 
PIM average, although slightly lower than the mean reported HERA-B cross section in Figs. |71 
andlHlfor the same values of m and /x, is still in good agreement. 
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Figure 9: The i^-factors for bb production in pp collisions with m = 4.75 GeV. We present Kq^^ 
(sohd), -ft'*^^-' (dashed) and -ft'sub (dot-dashed) results for (a) and (b), n = m/2, (c) and (d), 
^ = m and, (e) and (f), = 2m. Results with the GRV98 HO proton densities are shown on 
the left-hand side while results with the MRST2002 NNLO proton densities are shown on the 
right-hand side. 



The NNLO-NNNLL-t-C bb cross section at VS = 41.6 GeV with the MRST2002 NNLO 
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parton densities is 



_MRST2002NNLO _9oiQ+-'-^r,K 

"^NNLO-NNNLL+C — =■= J UU . I^O.iJ 

The central value represents the result for m = /i = 4.75 GeV, the first uncertainty is due to 
the scale variation and the second is due to the variation in m. The corresponding cross section 
for the GRV98 densities is somewhat smaller than the MRST result, 

GRV98 _ +^ , , ^ 

The uncertainties are also reduced for this density due to the lower value of A4 associated with 
this set. The NLO cross sections are, on the other hand, essentially identical for the two sets, 

+12 +10 , 

ctnlo = 17 _^ _g nb . (3.3) 

The K factors are shown for both sets of parton densities in Fig. El Recall that when the 
MRST2002 NNLO densities are used, all K factors are computed with these densities, as are 
the cross sections at each order shown in Fig. |H1 All the K factors computed with both densities 
tend to be larger for pp than ir'p collisions, especially at lower a/S". They are also stronger 
functions of y/s than those shown in Fig. El The value of Kq'^ is larger than K^"^^ except for 
fi = m/2. Note that -K'sub is less than K^^^ and fr*-^-* for all scales. We see that K^'^'' is nearly 
scale independent but has a value of ~ 2, still rather large. Including subleading terms gives 
K^^^ ^ 1 for > 20 GeV and /i = m/2, increasing to ~ 1.6 for = 2m. The values of K are 
nearly independent of the parton densities and bottom quark mass in all cases. 

Figure nm shows the scale dependence of bb production in pp collisions for both sets of parton 
densities. We again find a reduction of the scale dependence for the IPI NNLO-NNNLL+(^ 
results over all energies. Improvement in the approach of the ratio = m)/a{^ = /iq) toward 
unity is seen relative to the IPI NNLO-NNLL cross section ratio when /xq = m/2. The latter 
ratio is smaller than the NLO ratio, indicating stronger scale dependence for the IPI NNLO- 
NNLL cross section than that of the exact NLO result. Clear improvement is also seen for 
fiQ = 2m. 

In Fig. ^2 we plot the scale dependence for 0.3 < /x/m < 10 with = 41.6 GeV and 
m = 4.75 GeV. We show results for the Born, NLO, and NNLO-NNNLL+(^ cross sections. The 
scale dependence decreases with increasing order of the cross section. There is no plateau in 
the Born cross section while the exact NLO and the NNLO-NNNLL+C cross sections do show 
a peak in fi/m. The plateau at fi/m ^ 0.4 is broader for the NNLO-NNNLL+(^ cross section 
and the overall scale dependence is reduced relative to the exact NLO cross section. 

Figure ^1 shows the 6-quark transverse momentum distribution at HERA-B with y/S = 41.6 
GeV and m = 4.75 GeV. The Born, NLO, and NNLO-NNNLL+C results are shown on the left- 
hand side. On the right-hand side we plot K^^^ and K^"^^. We provide results with two different 



scales, fi = m and /i = m^ = yp^ + m^. There is some difference in the results with the two 
scales at larger values of pt- The distributions with the fixed scale, fi = m, are somewhat higher 
than those with /x = m^. Increasing pt and thus m^ decreases as, reducing da/dpr at higher 
Pt relative to the fixed scale choice. Although evolution increases the parton densities at low 
X, at the higher momentum fractions relevant here, evolution decreases the densities at higher 
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Figure 10: The scale dependence of bb production in pp collisions with m = 4.75 GeV. We give 
the ratios a{i^ = m)/a{n = fio) for the LO (dotted), NLO (solid), IPI NNLO-NNLL (dashed) 
and IPI NNLO-NNNLL+C (dot-dashed) cross sections. The results with //q = m/2 are given 
in (a) and (b) while the results with //q = 2m are given in (c) and (d). Results with the GRV98 
HO proton distributions are shown on the left-hand side while results with the MRST2002 
NNLO proton distributions are shown on the right-hand side. 
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Figure 11: The bottom quark cross section as a function oi fi/m in pp collisions with = 41.6 
GeV, m = 4.75 GeV and the MRST2002 NNLO densities. The Born (sohd), NLO (long- 
dashed), and NNLO-NNNLL+C (dashed) results are shown. 
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scales. Both the running of with scale and the evolution of the parton densities work to 
decrease the cross sections with the running scale relative to the fixed scale at higher p^- With 
both choices, we see an enhancement of the bottom quark transverse momentum distribution 
similar to that of the total cross sections when the threshold corrections are added. The shapes 
of the distributions, however, are very similar [Tl I25j. 

The K factors are, on average, similar to those shown in Fig. El at the same energy. When 
/i = m, both K factors tend to decrease with px over the range shown. As pr increases, the 
parton momentum fractions, x, probed also increases reducing both the available phase space 
and the gg contribution to the cross section. These two effects work together to reduce the K 
factors at larger pt [2^1 • When n = rriT, both K factors are essentially constant, as also seen 
in Ref. 




Pt (GeV) p, (GeV) 

Figure 12: The bottom quark transverse momentum distributions at = 41.6 GeV, m = 4.75 
GeV and the MRST2002 NNLO parton densities. On the left-hand side we show the Born 
(solid), NLO (long-dashed), and NNLO-NNNLL+(^ (dashed) results while on the right-hand 
side we present Kq^^ (solid) and K^'^^ (dashed). The bold lines are calculated with /i = rriT 
while the thin lines are with fi = m. 



4 Charm quark production 

We now turn to charm quark production. There is much more data on charm than bottom 
production. However, the early data are all rather low statistics and have not been updated 
to include the latest branching ratios. In our previous report on charm production to NNLO- 
NNLL we incorporated all the data as published, making assumptions only on how to 
extrapolate from the measured D meson cross sections to the charm quark cross sections. 
In this paper we choose only to compare our results to the n'p and pp data tabulated in 
Refs. fI7\ I2H] ■ These data are the most recent and incorporate the newest measurements of 
branching ratios. The 7r~p data in the threshold region, y/S < 30 GeV, in order of increasing 
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are found in Refs. 123 EOl ISH 1^21 while the pp data in the same energy region are given 
in Refs. j2niEniEll- 

How the cc pairs hadronize is a particularly important question for energies near threshold 
where some channels may be energetically disfavored. We follow Ref. and assume that 
a{Ds)/(y{D^ + D+) ^ 0.2 and a(Ac)/(T(D° + D+) ^ 0.3, independent of energy, so that the 
total cc cross section is obtained from ^ l.ba{DD). This assumption could have a strong 
energy dependence near threshold. Thus as many charm hadrons (mesons and baryons) as 
possible should be measured to better understand fragmentation and hadronization. Finally, 
some of the data are taken on nuclear targets and then extrapolated to 7r~p and pp assuming 
a linear A dependence jHSHSEI- 

Recent comparisons of the full pp — cc data set with exact NLO cross sections were made 
to determine the best mass and scale choices for extrapolation to higher energies [H7]. Rough 
agreement with the data up to the top ISR energy, = 63 GeV, was found for m = 1.2 GeV 
and IX = 2m for the MRST densities and m = 1.3 GeV with fi = m for the GRV98 densities 
j37j . These values of m are rather small compared to the typical value of 1.5 GeV. Thus, as in 
our previous paper we calculate the NLO, IPI NNLO-NNLL and, in addition here, the IPI 
NNLO-NNNLL+(^ cross sections using m = 1.2, 1.5 and 1.8 GeV as well as fi = m and 2m. 
We can then test whether the NNLO+NNNLL+(^ cross sections might favor a higher charm 
quark mass. Our charm calculations also employ the GRV98 HO and MRST2002 NNLO proton 
parton densities. 

4.1 Charm production in tt j9 interactions 

We first discuss the results from vr^p interactions. These data are typically reported as the 
forward cross sections, xp > 0. The cross section is extrapolated to all xp assuming that 
o"(all xp) = l.Qcr{xF > 0), calculated to NLO. The difference is not a factor of two as in pp 
production since the forward xp distribution is harder than the backward distribution due to 
the harder large x behavior of the pion gluon distribution. 

The reaction 7i~p — >■ cc is dominated by the gg channel for > 15 GeV ^H]- The gg and 
qq channels are only similar in magnitude very close to threshold. The comparison of the data 
with the exact NLO, IPI NNLO-NNLL and IPI NNLO-NNNLL+C cross sections, calculated 
with the GRV98 HO proton parton densities and the GRS pion parton densities, is shown in 
Fig. El The exact NLO and the IPI NNLO-NNLL results were previously shown in Ref. j^. 
The new IPI NNLO-NNNLL+C results, indicated by the dot-dashed curves, are similar to the 
average of the NNLO-NNLL IPI and PIM cross sections of Ref. |3]. At the higher end of the 
"s/S* range studied, the new results with the subleading logs are higher than this average because 
the PIM NNLO-NNLL gg contribution is large and negative above threshold. Eventually this 
contribution becomes larger than the exact NLO cross section, resulting in a negative PIM 
NNLO-NNLL total cross section and a decreased IPI and PIM average NNLO-NNLL cross 
section, particularly for the lower values of m. Inclusion of the subleading logs mitigates 
this behavior [S] albeit not enough to make the PIM calculation reliable for gluon dominated 
processes. 

There is relatively good agreement between the exact NLO calculations and the data for 
m = 1.2 GeV and fi = 2m, Fig. IT^ b). We can see, however, that the data are also relatively 
compatible with the IPI NNLO cross sections, both the NNLO-NNLL and NNLO-NNNLL+C 
calculations, when m = /i = 1.5 GeV. This result suggests that a full NNLO calculation would 
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Figure 13: The energy dependence of cc production in 7r~p collisions with, (a) and {h),m — 1.2, 
(c) and (d), m = 1.5, and, (c) and (f), m = 1.8 GeV, calculated with the GRV98 HO proton 
densities and the GRS pion densities. We show the NLO (solid), IPl NNLO-NNLL (dashed) 
and IPI NNLO-NNNLL+^ (dot-dashed) results. On the left-hand side, /j, = m while on the 
right-hand side, /i — 2m. 
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be more compatible with a larger charm quark mass. 
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Figure 14: The fT-factors for cc production in tt p collisions with, (a) and (b), m = 1.2, (c) 
and (d), m = 1.5, and, (e) and (f), m = 1.8 GeV. We present K^^^ (sohd), K^"^^ (dashed) and 

(2) 

K^J^ (dot-dashed) for = m (left-hand side) and fi = 2m (right-hand side). 



The 7r~p charm production K factors are shown in Fig. As expected, the K factors are 
all larger for charm than bottom production due to the dominance of the gg channel over all 
as well as the larger value of a^. There is a significant decrease in the K factor when the 
subleading terms are included: K^'^^ ^ 1.5 whereas K'^'^^ ^ 2. All the K factors decrease slightly 
with increasing m. The mass effect is smaller on K^^^ than on K^^\ the IPI NNLO-NNLL K 
factor. Indeed the K^'^^ mass dependence is larger than either that of Kq^'^ or K^'^^^. The exact 
NLO K factor, Kq^\ shows the strongest scale dependence, increasing from ^ 2 for fi = m to 
f» 2.5 when n = 2m, as also seen for bb production. Note also that while all the K factors are 
relatively energy independent, K^^^^ exhibits the largest a/S* dependence. 

Figure 1131 demonstrates the scale dependence of charm production in 7r~p interactions. Since 
/X = m/2 is less than the initial scale of most sets of parton densities, this value of is essentially 
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Figure 15: The scale dependence of cc production in tt p collisions with (a) m = 1.2, (b) m = 
1.5, and (c) m = 1.8 GeV. We give the ratios a{fi = m)/a{fi = 2m) for the LO (dotted), NLO 
(sohd), IPI NNLO-NNLL (dashed) and IPI NNLO-NNNLL+C (dot-dashed) cross sections. 
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meaningless for charm production. Therefore, we plot the ratio cr(/i = m)/a{fi = 2m) as a 
function of for all three values of m. The scale dependence only weakly depends on charm 
mass although the dependence is slightly stronger at low v^- The IPI NNLO-NNNLL+(^ ratio 
shows significant improvement over both the NLO and the IPI NNLO-NNLL results. 



4.2 Charm production in pp interactions 




Figure 16: The energy dependence of cc production in pp collisions with, (a) and (b), m = 1.2, 
(c) and (d), m = 1.5, and, (e) and (f), m = 1.8 GeV, calculated with the GRV98 HO proton 
densities. We show the NLO (sohd), IPI NNLO-NNLL (dashed) and IPI NNLO-NNNLL+C 
(dot-dashed) results. On the left-hand side, fi = m while on the right-hand side, fi = 2m. 

We now consider pp — > cc interactions. In Figs. El and El "we compare the exact NLO, IPI 
NNLO-NNLL and IPI NNLO-NNNLL+C cross sections calculated with the GRV98 HO and 
MRST2002 NNLO proton parton densities respectively. There is a larger difference between 
the results with the two parton densities than seen for bb production. Since A3 > A4 the overall 
NNLO corrections are thus larger for charm than bottom as well, as evident from the larger 
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charm K factors. At NLO, the best agreement is again with m = 1.2 GeV and /i = 2m, seen 
in Figs. ITHl and IT7I (b). The new result in Fig. ^Jis the subleading NNLO-NNNLL+(^ cross 
section, the other calculations, also shown in Ref. jlj, are repeated for comparison purposes. 
The results in Fig. El are all new. 

The MRST2002 NNLO parton densities generally give larger cross sections, even for the 
exact NLO result since the value of A3 is larger than that of the GRV98 HO set. Also due the 
larger A3, the NNLO corrections are significantly larger, as seen in Fig. El A similar effect was 
seen in the comparison of the GRV98 HO and CTEQ5M results in Ref. j^. 




VS [GeV] VS [GeV] 

Figure 17: The energy dependence of cc production in pp collisions with, (a) and (b), m = 1.2, 
(c) and (d), m = 1.5, and, (e) and (f), m = 1.8 GeV, calculated with the MRST2002 NNLO 
proton densities. We show the NLO (sohd), IPI NNLO-NNLL (dashed) and IPI NNLO- 
NNNLL+C (dot-dashed) results. On the left-hand side, fi = m while on the right-hand side, 
/i = 2m. 

In both Figs. ITHl and IT71 the exact NLO cross sections calculated with m = 1.2 GeV and 
/i = 2m are relatively compatible with the data. Note that in Fig. ITCT b). the IPI NNLO- 
NNNLL-|-^ result is in somewhat better agreement with the two highest energy data points 
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than the exact NLO. The IPI NNLO cross sections are in rather good agreement with the 
data when m = fi = 1.5 GeV is used with the MRST2002 NNLO parton densities. Indeed, 
the IPI NNLO-NNNLL+(^ result in Fig. ITTT c) agrees rather well with the two higher energy 
data points. Thus, as in Ref. 0, we can conclude that the full NNLO result can likely describe 
the charm data well with m = = 1.5 GeV whereas the lower mass is needed with an NLO 
calculation alone. 
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Figure 18: The i^'-factors for cc production in pp collisions with, (a) and (b), m = 1.2, (c) and 
(d), m = 1.5, and, (e) and (f), m = 1.8 GeV, calculated with the GRV98 HO parton densities. 
We present Kq^^ (solid), K^'^'^ (dashed) and k!:^^ (dot-dashed) for /i = m (left-hand side) and 
fi = 2m (right-hand side). 

Table ^ gives the charm cross sections in pp collisions at two energies, ^/S = 6.98 GeV, 
the future GSI energy, and 17.3 GeV, the CERN SPS ion energy. We present results for both 
the MRST2002 NNLO and GRV98 parton densities based on a central value of m = yU = 1.5 
GeV. (The choice of mass and scale used for our central value is for better illustration of the 
uncertainties rather than any fit to data.) The first uncertainty is due the the scale choice. 
Since we do not calculate the result for /z = m/2 here, we show only the difference between the 
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values of /i = m and 2m. The second uncertainty is that due to the charm quark mass. The 
exact NLO and the IPI NNLO-NNNLL+(^ cross sections are shown. Note that the differences 
between the MRST2002 and GRV98 densities are larger for charm since the difference in the 
values of A3 between the two sets is larger than that for A4. The relative increase of the NNLO- 
NNNLL+C cross section is highest at the lower energy with an increase over the NLO result 
by a factor of 2.6 for the MRST densities and 2.2 for the GRV98. At 17.3 GeV, the increase is 
by a factor of 1.8 and 1.5 respectively. Thus the importance of the near-threshold corrections 
is reduced further from the charm production threshold, as might be expected. 



a (/ib) 


(GeV) 


Order 


MRST2002 NNLO 


GRV98 


6.98 


NLO 


_l_n c;r 

0.034 -0.027 ^(, (,32 


-i-n 42 

0.028 - 0.022 

—0.026 


6.98 


NNLO-NNNLL+C 


0.09 - 0.07 «4 


0.061 - 0.05 


17.3 


NLO 


3.8 -2.1 


2.8 -1.4 


17.3 


NNLO-NNNLL+C 


fi7 '^A +22-5 
6.7 -3.4 _^_g 


4.1 -1.8 



Table 1: The cc production cross sections in pp collisions at \^ = 6.98 and 17.3 GeV. The 
exact NLO results and the approximate NNLO-NNNLL+C results, based on m = /i = 1.5 GeV, 
are shown. The first uncertainty is due to the scale choice, the second, the charm quark mass. 

In Figs. CHI and [H we compare the pp K factors for GRV98 HO and MRST2002 NNLO 
respectively. The same trends are seen as in tt~p interactions although the energy dependence 
at low -s/S, particularly for K'^^\ is stronger. None of the K factors are strong functions of 
mass, scale or parton density. The parton density dependence is thus reduced relative to the 
calculations with CTEQ5M shown in Ref. This is perhaps because the CTEQ5M gluon 
density is higher than either the GRV98 HO or MRST2002 NNLO gluon densities at x < 0.1 and 
at the low charm production scales. Although the initial scales of the CTEQ5 and MRST sets 
are not very different, 1 GeV^ for CTEQ5M and 1.25 GeV^ for MRST2002, the low x behavior is 
quite different. Indeed at low scales, when m = = 1.2 GeV, the CTEQ5M gluon distribution 
is essentially constant as x — > with the tail of the distribution appearing at x ~ 0.1, while 
for the MRST2002 NNLO set the gluon distribution becomes negative at x < 0.001 and peaks 
at X < 0.1. When /i is increased to 2m for the same mass, the CTEQ5M density decreases 
with increasing x while the MRST2002 NNLO density peaks around x ~ 0.001, decreasing at 
lower X. Although our calculations are not sensitive to the very low x region, the different 
low-x behavior affects the gluon distributions at higher x. The lower MRST gluon density, 
combined with the lower A3 relative to CTEQ5M reduces the higher order corrections and 

('2) 

hence the K factors. Thus K^J^^ varies between 1.3 and 1.7 for GRV98 HO and 1.5 to 1.9 for 
MRST2002 NNLO when > 15 GeV, significant improvements over the IPI NNLO-NNLL 
resuhs, -2-2.4 for GRV98 HO and 2.5 - 3 for MRST2002 NNLO and CTEQ5M. In 
addition, K'^y^ < Kq^^ for all cases considered, not a feature of the IPI NNLO-NNLL result. 
We compare the scale dependence of the cross sections in Fig. 1201 The GRV98 HO cross 
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Figure 19: The X-factors for cc production in pp collisions with, (a) and (b), m — 1.2, (c) and 
(d), m = 1.5, and, (e) and (f), m = 1.8 GeV, calculated with the MRST2002 NNLO parton 
densities. We present Kq^'' (solid), K^^^ (dashed) and K^^^ (dot-dashed) for = m (left-hand 
side) and /i — 2m (right-hand side). 
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Figure 20: The scale dependence of cc production in pp collisions with, (a) and (b), m = 1.2, 
(c) and (d), m = 1.5, and, (e) and (f), m = 1.8 GcV. We give the ratios a{fi = m)/a{fi = 2m) 
for the LO (dotted), NLO (solid), IPI NNLO-NNLL (dashed) and IPI NNLO-NNNLL+C (dot- 
dashed) cross sections. Results with the GRV98 HO densities are given on the left-hand side 
while the MRST2002 NNLO results are shown on the right-hand side. 
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section ratios on the left-hand side are compared to the MRST2002 NNLO ratios on the right- 
hand side. The scale dependence is similar for the two sets of parton densities although the 
MRST scale dependence is somewhat stronger at low v^- The IPI NNLO-NNNLL-t-^ scale 
dependence is reduced relative to the IPI NNLO-NNLL dependence which is stronger than the 
exact NLO dependence for most values of mass. Indeed the IPI NNLO-NNLL ratios are only 
similar to the NLO ratios for m = 1.8 GeV while the IPI NNLO-NNNLL+C ratios are lower 
than the rest of the calculated ratios for all masses. 

In Fig.|2I]we plot the scale dependence for 0.8 < /x/m < 10 at ^/S = 17.3 GeV and m = 1.5 
GeV. We show results for the Born, NLO, and NNLO-NNNLL+(^ cross sections. None of the 
cross sections show a plateau at any value of /i/m, as seen for the NLO and NNLO bb cross 
sections in Fig. [TT] However, the relative dependence of the cross section on /i/m decreases as 
the accuracy of the cross section increases. 




Figure 21: The charm production cross section as a function of /i/m at = 17.3 GeV with 
m = 1.5 GeV and the MRST2002 NNLO parton densities. We show the Born (solid), NLO 
(long-dashed) and NNLO-NNNLL+C (dashed) results. 

Figure |221 shows the charm quark transverse momentum distributions in pp collisions at 
= 17.3 GeV, m = 1.5 GeV, and two scales, fi = m and /i = = ^p"^ + m?. The Born, 
NLO, and NNLO-NNNLL+C results are shown on the left-hand side. On the right-hand side 
we plot k'^'^ and k'^^-^. a gain the enhancement in the pt distribution is similar to that seen for 
the total cross sections in Fig. The behavior of the K factors as a function of px is similar 
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to that in Fig. [H 

8 I 6 




Pt (GeV) p, (GeV) 

Figure 22: The charm quark transverse momentum distributions at = 17.3 GeV and 
m = 1.5 GeV with the MRST2002 NNLO parton densities. On the left-hand side we show the 
Born (sohd), NLO (long-dashed), and NNLO-NNNLL-|-^ (dashed) results while on the right- 
hand side we present K^^^ (solid) and (dashed). The bold lines are calculated with /z = rriT 
while the thin lines are with fi = m. 



5 Conclusions 

In this paper we have calculated soft NNLO corrections to the bottom and charm quark total 
cross sections and transverse momentum distributions in hadron-hadron collisions. We have 
added new subleading soft NNNLL terms and some virtual terms, including all soft-plus- virtual 
factorization and renormalization scale-dependent terms. We have found that these new sub- 
leading corrections reduce the size of the NNLO cross sections, and thus the K factors, as well 
as diminish the scale dependence of the cross section. 
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